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FOREWORD 


This  report  consists  of  the  second  phase  of  an  analytical  program 
to  provide  fundamental  solutions  to  the  non-uniform  flow  fields  of  wing- 
propeller  aerodynamics,,  This  program  was  originated  by  Dr,  Scott 
Rethorst  as  an  extension  of  his  basic  lifting  surface  solution  obtained  in 
an  earlier  study  conducted  at  the  California  Institute  of  Technology,  ^ 

The  program  was  formulated  so  that  the  effects  of  primary  variables 
such  as  wing  plan  form  were  first  determined.  The  analysis  is  then 
extended  in  order  of  increasing  complexity  to  enhance  its  generality  and 
to  encompass  secondary  features.  Each  phase  of  the  program  establishes 
the  foundation  for  the  succeeding  phases  in  an  orderly  growth  pattern. 

The  following  four  phases  constitute  the  basic  analytical  portion  of  this 
program: 

1,  Single  Jet  Theory  -  The  Rethorst  analysis  of  Ref,  1 
was  applied  to  an  extensive  systematic  spectrum  of  43  7  wings.  Optimum 
finite  wing  planforms  were  determined,  ^ 

2,  Multiple  Jet  Theory  -  The  single  jet  theory  is  further 
generalized  and  applied  to  the  case  of  multiple  jets.  The  effect  of  slip¬ 
stream  rotation  is  also  included,  (Present  Report), 

3,  High  Angle  of  Attack  Theory  -  The  generalized  theory 
is  applied  to  analyses  of  (1)  wings  located  at  various  heights  in  the  jet, 

(2)  highly  cambered  wings  as  used  in  deflected  slipstream  V/STOL  con¬ 
cepts,  and  (3)  tilt  wing  configurations  where  the  jet  is  at  an  angle  to  the 

(3) 

free  stream  flow, 

4,  Separated  Flow  Theory  -  Hydrodynamic  cavitation  theory 


is  used  to  determine  the  separated  flow  field  incurred  during  the  transi¬ 
tional  flight  of  V/STOL  aircraft.  ^ 

The  original  Rethorst  analysis  of  a  wing  extending  through  a  single 
jet  as  used  in  Phase  1  employs  an  infinitesimal  vortex  element,  which 
early  in  the  analysis  was  integrated  over  a  finite  span  element.  The 
wing  was  then  represented  by  a  series  of  such  finite  spanwise  horseshoe 
vortex  elements. 

Each  such  horseshoe  vortex  element  was  then  split  into  components 
which  were  even  and  odd  with  respect  to  the  direction  of  flow.  The  flow 
field  due  to  the  two  dimensional  even  component  was  found  in  a  straight¬ 
forward  algebraic  manner.  The  flow  field  due  to  the  three  dimensional 
odd  component  is  quite  complex,  involving  infinite  sums  of  infinite  integrals 
of  Bessel  functions.  Electronic  data  processing  equipment  is  necessary 
for  evaluation  of  these  quadratures. 

The  present  paper,  in  extending  the  Rethorst  analysis  to  multiple 
jets,  introduces  four  further  significant  refinements.  First,  the  infinitesi¬ 
mal  vortex  elements  are  retained  throughout  the  analysis,  delaying  the 
integration  over  the  span  until  a  final  step.  This  technique  simplifies 
the  analysis  and  improves  the  accuracy  of  the  results.  Second,  the  sym¬ 
metry  restrictions  have  been  removed  so  that  more  generalized  geometries 
may  be  analyzed.  Third,  the  odd  component  has  been  further  partitioned 
into  a  simple  algebraic  component  plus  a  complex  remainder  containing 
all  of  the  Bessel  function  terms.  If  this  remainder  can  be  bounded  and 
shown  to  be  small  enough  to  be  disregarded,  even  in  certain  special  cases, 
a  major  saving  in  computation  will  be  achieved.  Fourth,  the  effects  of 
slipstream  rotation  have  been  included  as  a  super-imposed  element. 


SUMMARY 


The  basic  Re'.horst  lifting  surface  solution  of  a  wing  extending 
through  a  single  jet  was  generalized  to  enhance  its  applicability  to  the 
solution  of  many  general  and  secondary  problems  concerned  with  the 
non-uniform  aerodynamics  of  wing-propeller  interaction.  The  general¬ 
ization  is  first  developed  in  detail  in  terms  of  single  jet  theory.  Then 
it  is  applied  to  multi-jet  arrangements. 

Several  significant  refinements  to  the  original  Rethorst  theory 
are  introduced.  The  wing  is  represented  by  a  distribution  of  infinitesi¬ 
mal  vortex  elements  instead  of  a  large  but  finite  number  of  horseshoe 
vortices.  Certain  symmetry  restrictions  are  removed  and  the  effect 
of  slipstream  rotation  is  included. 

The  spanwise  lift  distribution  of  an  example  multi-jet  arrange¬ 
ment  was  determined.  Large  increases  in  lift  inboard  and  within  the 
jets  were  obtained  when  the  free  stream  velocity  was  small  compared 
to  the  jet  velocity.  A  lesser  but  significant  increase  occurred  outboard 
of  the  jets. 

This  large  lift  magnification  offers  the  potential  of  large  im¬ 
provements  in  STOL  capabilities  and  warrants  a  further  computational 
effort  which  is  presently  being  conducted. 
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I.  INTRODUCTION 


The  problem  of  the  spanwise  lift  distribution  on  a  wing  of  finite 

aspect-ratio  extending  through  a  central  jet  slipstream  has  been  investi- 

(5) 

gated  by  several  authors.  Koning  first  treated  this  problem,  using  the 
lifting  line  theory  together  with  certain  simplifying  assumptions.  In 
Koning's  work,  the  slipstream  is  taken  to  be  a  circular  jet  of  air  moving 
with  a  uniform  velocity  parallel  to  the  main  stream.  The  distortion  of 
the  slipstream  due  to  the  presence  of  the  wing  is  assumed  to  be  small 
so  that  the  boundary  conditions  on  the  jet  boundary  may  be  applied,  by 
linearization,  at  its  undisturbed  position.  Furthermore,  the  difference 
between  the  jet  velocity  and  the  main  flow  velocity  is  assumed  small; 
this  limitation  was  introduced  primarily  to  simplify  the  method  of 
calculation  used  by  Koning.  Recently  Koning's  theory  has  been  extended 
to  render  results  valid  for  all  jet  velocities  (with  the  other  original 
assumptions  retained)  by  Graham,  Lagerstrom,  Licher  and  Beanei^, 

The  lift  distribution  along  a  wing  of  finite  span  with  a  central 

(7) 

jet  has  been  calculated  by  Squire  and  Chester  ,  using  the  lifting  line 
theory  together  with  a  rather  elaborate  expansion  of  the  circulation 
distribution  along  the  wing0  Although  the  lifting  line  approximation 
well  represents  a  wing  of  high  aspect-ratio  in  a  uniform  main  flow, 
the  same  approximation  adapted  to  the  present  problem  will  further 
require  that  the  slipstream  diameter  be  large  compared  with  the  wing 
chord.  In  practice,  however,  the  slipstream  diameter  and  the  chord 
are  usually  of  the  same  order  in  magnitude.  The  applicability  of  the 
lifting  line  theory  is  therefore  not  well  established. 
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On  the  other  hand,  if  the  slipstream  diameter  is  small  compared 
with  the  wing  chord,  the  lift  increase  due  to  the  jet  effect  is  produced 
mainly  over  a  low  aspect-ratio  surface,  namely  the  portion  of  the  wing 
in  the  jet  and  its  immediate  neighborhood.  This  special  case  has  been 
treated  by  applying  the  slender  body  theory  in  Ref,  6.  In  this  theory  the 
basic  velocity  potential  <pg  of  the  wing  alone  without  the  jet  is  considered 
as  known  and  the  jet  perturbation  potential  A <f>  =  <f>  -  <j>  is  assumed  to 
satisfy  the  slender  body  approximations;  the  problem  of  Atj>  is  then 
solved  with  the  effect  of  finite  aspect-ratio  on  A<fr  neglected. 

Judging  from  the  different  assumptions  inherent  in  these  two 
theories,  it  is  expected  that  the  lifting  line  theory  will  be  valid  for  large 
values  of  the  ratio  of  jet  diameter  to  wing  chord,  whereas  the  slender 
body  theory  will  be  valid  for  small  values  of  this  ratio.  In  order  to  as¬ 
certain  the  range  of  validity  of  these  two  theories,  their  results  have 
been  compared  in  Ref.  6  for  a  particular  example  (a  wing  with  sinusoidally 
varying  incidence  along  the  span,  extending  through  a  plane-sided  free 

jet  of  finite  width  and  infinite  height)  with  the  result  of  a  third  method, 

(81 

namely  the  approximate  lifting  surface  theory  of  Weissinger'  . 

The  comparison  shows  that,  at  least  for  this  particular  case,  the 
slender  body  theory  is  in  good  agreement  with  the  lifting  surface  theory 
when  the  jet  width  is  less  than  the  wing  chord.  The  lifting  line  theory 
is  inadequate  in  this  range  of  width-to-chord  ratio,  but  in  the  limit  tends 
to  the  Weissinger  theory  for  large  values  of  this  ratio  (about  10  or  larger). 
While  this  finding  indicates  the  range  of  validity  of  the  two  limiting  theories 
in  this  special  case,  their  respective  ranges  of  validity  in  the  general 
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case,  however,  are  not  definitely  established,  as  indicated  by  the  com¬ 
parison  (given  below)  of  their  results  with  the  existing  experimental  data. 

For  the  general  case  it  seems  evident  that  some  kind  of  lifting 
surface  theory  would  be  required  to  describe  properly  the  actual  flow. 
Based  on  this  point  of  view,  Rethorst^  treated  the  original  Koning 
problem  by  using  the  Weissinger  approximation.  In  this  analysis,  the 
wing  is  represented  by  a  large  number  of  horseshoe  vortices,  which 
are  obtained  by  integration  of  an  infinitesimal  vortex  over  a  small  but 
finite  span  element.  The  problem  is  then  solved  by  the  finite  step  method. 
To  avoid  the  difficulties  encountered  in  evaluating  the  convergence  of  the 
solution  at  the  jet  boundary,  no  vortex  is  chosen  to  stride  across  the 
jet  boundary.  This  arrangement  simplifies  the  boundary  condition  on 
the  downwash  at  the  intersection  of  the  jet  boundary  and  the  wing.  The 
lift  distribution  given  by  Rethorst's  theory  is  shown  to  lie  between  the 
two  limiting  theories^.  This  theory  of  Rethorst  has  been  subsequently 
applied  to  compute  various  flow  quantities  for  a  wide  range  of  the  velocity 
ratio  (main  stream  velocity  to  jet  velocity)  and  the  ratio  of  jet  diameter 
to  chord;  the  results  are  reported  in  Ref,  2. 

There  are  only  a  few  experimental  results  on  this  slipstream 

(9) 

problem.  A  series  of  experiments  were  conducted  by  Smelt  and  Davies  , 
using  a  rectangular  wing  held  at  an  incidence  of  6.  6  degrees,  and  with 
the  propeller  axis  parallel  to  the  main  flow.  Only  the  aerodynamic  forces 
on  the  wing  itself  were  measured.  Compared  with  this  set  of  experiments, 
the  lift  increment  predicted  by  the  slender  body  theory  is  about  30%  lower 
than  the  experimental  value^,  whereas  the  lift  increment  given  by  the 
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(1) 


lifting  line  theory  is  more  than  double  the  experimental  data. 

In  a  set  of  experiments  carried  out  by  Stuper^^,  the  spanwise 
lift  distribution  on  a  rectangular  wing  was  measured  at  two  angles  of 
attack.  A  specially  designed  propeller  was  used  to  produce  a  jet  with 
nearly  uniform  velocity  distribution  and  without  rotation.  A  comparison 
shows  that  the  experimental  lift  distribution  of  Stuper  lies  between  the 
two  limiting  theories,  whereas  Rethorst's  theory  is  found  in  substantial 
agreement  with  the  experiment  within  the  jet.  Outside  the  jet,  the 
correlation  is  obscured  by  the  use  of  an  upstream  duct  in  Stuper's  experi¬ 
ment. 

Another  experimental  study  of  this  problem  was  made  recently 
by  Brenckman^^.  It  was  found  that  within  a  range  of  moderate  wing 
incidences,  a  substantial  part  of  the  lift  increment  produced  by  the  slip¬ 
stream  is  attributed  to  a  destalling  or  boundary-layer-control  effect. 

With  this  part  of  the  lift  increment  deducted,  the  remaining  part  was 
founa  in  good  agreement  with  the  slender  body  theory.  The  slipstream 
velocity  profile  (deviation  from  a  uniform  distribution),  rotation,  and 
deformation  (from  its  original  circular  shape)  were  claimed  to  be  of 
secondary  importance.  However,  no  apparent  attempt  was  mentioned 
in  Brenckman's  investigation  to  examine  the  case  when  the  angle  of 
attack  is  so  small  that  the  destalling  effect  is  practically  absent,  and 
to  examine  the  other  ranges  of  the  chord-to-diameter  ratio.  Thus  it 
seems  that  further  experimental  investigation  of  these  points  would  be 
required  in  order  to  provide  a  more  sound  basis  for  comparison  between 
the  different  theories  and  experiments. 
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More  recently  an  extensive  experimental  investigation  of  the 

slipstream  problem  was  undertaken  at  Princeton  University,  Department 

of  Aeronautical  Engineering.  The  first  report  of  this  work  is  on  finite 

(12) 

spans  and  produced  experimental  data  that  was  in  close  agreement 
with  the  Rethorst  theory, 

A  further  recent  Princeton  report  on  this  program  dealt  with 

(13) 

infinite  aspect  ratio  wings  .  The  Rethorst  formulation  again  proved 
to  yield  the  best  correlation  with  experimental  results. 

This  Princeton  program  provides  substantial  experimental  evi¬ 
dence  of  the  validity  of  the  present  theoretical  program,  particularly 
the  first  phase  based  on  Rethorst's  work.  This  favorable  correlation 
lends  encouragement  to  the  present  generalization  of  this  initial  effort. 

The  present  study  is  concerned  mainly  with  the  generalization 
of  the  Rethorst  lifting  surface  theory  so  that  the  analysis  may  be  more 
readily  applied  to  an  extensive  set  of  general  problems,  such  as  multiple 
jets,  slipstream  rotation,  ground  effect,  and  the  optimum  wing  form 
for  total  lift.  The  first  two  problems  are  solved  in  the  present  paper. 

Because  the  Weissinger  method  is  less  time-consuming  than  other 
more  elaborate  lifting  surface  theories,  it  is  adopted  here  as  a  first  ap¬ 
proximation  to  the  lifting  surface  theory.  In  order  to  simplify  the  analysis 
and  to  improve  the  accuracy  of  the  results,  the  wing  is  represented  by 
a  piece-wise  continuous  distribution  of  circulation  T,  which  may  admit 
a  step- jump  only  at  the  jet  boundaries.  The  circulation  distribution  T 
is  further  expressed  by  different  Fourier  series  in  different  regions 
(inside  and  outside  the  jets)  over  the  wing.  With  this  expansion  of  T  the 
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condition  that  the  downwash  at  the  jet  boundaries  must  be  bounded  can 
easily  be  enforced.  The  first  few  Fourier  coefficients  are  determined 
by  applying  a  number  of  boundary  conditions  at  approximate  points  on 
the  wing. 

The  case  of  a  single  jet  is  first  considered  in  detail;  the  analysis 
is  then  generalized  to  the  case  of  multiple  jets.  The  effect  of  the  slip¬ 
stream  rotation  is  accounted  for  in  an  approximate  manner. 
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II.  ANALYSIS 


A.  GENERAL  FORMULATION 


In  this  analytical  treatment  of  the  effect  of  a  number  of  slip¬ 
streams  on  the  spanwise  lift  distribution  over  a  wing,  the  following 
simplifying  assumptions  are  introduced.  Each  slipstveam  is  taken  to 
be  a  jet  of  air  which  has  a  circular  cross  section  and  has  its  central 
axis  parallel  to  the  main  stream  surrounding  the  jets.  To  begin  with, 
the  slipstreams  are  assumed  to  have  no  rotation  about  their  axes;  the 
effect  of  small  slipstream  rotation  will  be  approximately  estimated 
later  in  the  analysis.  Thus  the  slipstream  velocity  and  the  main  stream 
velocity  are  assumed  to  be  originally  uniform  and  constant,  equal  to  Vj 

and  V  respectively.  In  this  analysis  the  quantity  (V.- V  )/(V.+  V  ) 

J  J 

need  not  be  limited  to  be  small.  When  V.  »  Y  ,  the  effect  of  the  slip- 

J  o 

streams  becomes  predominant;  consequently  the  assumption  that  the 
slipstreams  are  parallel  to  the  main  stream  becomes  less  important 
and  may  thereby  be  relaxed  to  some  extent.  It  may  also  be  noted  that 
for  Vj  »  Vq,  the  mixing  region  at  the  jet  boundary  due  to  the  viscous 
effect  will  become  appreciable;  this  real  fluid  effect,  however,  will  be 
neglected.  It  should  be  mentioned  that,  as  pointed  out  by  Brenckman^^, 


a  substantial  part  of  the  lift  increment  at  moderate  incident  angles  pro¬ 
duced  by  the  slipstream  is  attributed  to  a  "destalling"  effect  which  reduces 
the  flow  separation  on  account  of  the  high  velocity  slipstream.  This  point 
of  view  would  seem  to  require  some  kind  of  new  or  modified  analysis, 
presumably  with  the  use  of  the  free  streamline  theory  to  represent  the 
separated  flow.  Within  the  present  framework,  however,  the  flow  is 
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assumed  to  be  free  from  separation 

The  flow  is  assumed  to  be  incompressible  and  nonviscous,  hence  it  pos¬ 
sesses  a  velocity  potential  inside  and  outside  the  jet,  though  the  potential 
need  not  be  continuous  at  the  jet  boundary.  The  jet  boundary  may  be  regarded 
as  a  vortex  sheet  separating  the  inside  and  outside  flow.  The  only  other 
singularities  are  the  vortex  elements  representing  the  wing  and  the  trail¬ 
ing  vortex  sheet  shed  from  the  wing.  The  disturbances  produced  by  the 
wing  are  assumed  to  fall  off  far  away  from  the  wing  except,  of  course, 
along  its  trailing  vortex  sheet. 

Before  the  analysis  is  extended  to  the  case  of  multiple  jets,  we 
consider  first  the  case  of  a  single  circular  jet,  with  its  radius  normalized 
to  unity  for  simplicity.  The  jet  axis  is  taken  to  be  the  x-axis,  directed 
along  the  main  flow  and  centered  at  the  wing  chord.  The  wing  extends 
through  the  jet,  spanning  from  y  =  -b^  to  y  =  +b^  (b^,  b^  >  1,  see  Fig.  1). 

The  z-axis  is  taken  to  point  upward.  The  ratio  between  the  main  stream 
velocity  Vq  and  the  jet  velocity  Vh  will  be  denoted  by 

p  =  v0/v..  0  <  p  <  1  .  (1) 


The  perturbation  velocity  potential  outside  and  inside  the  jet  will  be  de¬ 
noted  by  cp ^  and  cpj  so  that  the  total  velocity  will  be  respectively 


%  =  VT  +  grad  (pQ,  q.  =  vT  +  grad  cp. 


(2) 


where  i  is  a  unit  vector  along  the  x-axis.  In  their  respective  flow 

regions  cp  and  (p.  satisfy  the  Laplace  equation 
J 

VZ(P0  =  o  ,  vVj  =  o  .  (3) 


The  distortion  of  the  slipstream  by  the  wing  is  assumed  small  so  that 
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by  linearization,  the  conditions  on  the  jet  boundary  may  be  applied  on 

the  original  circular  free  boundary.  There  are  two  boundary  conditions 

at  the  jet  boundary.  One  of  them  is  dynamic  in  nature,  stating  that  the 

(5) 

pressure  must  be  continuous  at  the  jet  boundary;  that  is, '  ' 

<pj  =  V-(p0  on  r  =  (y2+  z2)1/2  =  1  .  (4) 

The  other  condition,  kinematic  in  nature,  expresses  that  the  inside  and 
outside  flow  must  be  tangential  to  each  other  at  the  boundary,  that  is, 

9<P0/9r  =  H9<Pj/9r  on  r  =  1  (5) 

The  boundary  conditions  on  the  trailing  vortex  sheet  of  the  wing  will  be 
described  below  (immediately  after  Eq.  6),  and  the  boundary  condition 
on  the  lifting  surface  representing  the  wing  will  be  given  later  in  this 
analysis  (see  Section  B4), 

In  order  to  generalize  the  analysis  and  to  improve  its  accuracy 
for  subsequent  numerical  computations,  the  wing  will  be  represented  by 
a  distribution  of  infinitesimal  vortex  elements  instead  of  a  number  of 
finite  horseshoe  vortices  as  in  Ref,  1,  Let  us  first  consider  an  infinitesi¬ 
mal  element  of  bound  vortex  of  strength  I"’('n)6'n  centered  at  x  =  0,  y  =  r|, 
z  =  0,  and  extending  over  a  line  element  6r|  along  the  y-axis,  the  free 
stream  being  everywhere  uniform,  of  velocity  V  in  the  positive  x-direc- 
tion.  The  trailing  vortex  sheet  of  this  bound  vortex  element  lies  approxi¬ 
mately  in  the  region  x  >  0,  r|  -  6r(/2  <  y  <  r|  +  6t)/2,  z  =  0.  The  solution 
of  this  problem  is  well  known  (e,  g,  see  Ref.  14);  it  may  be  expressed  as 

5 (p  =  5<^1+  6 (pz  =  ^  r(ri)6r)[  F^y-t],  z)  +  F2(x,  y-t],  z)j,  (6a) 
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where 


Fl(y*  z)  =  ~2 - 7 

y  +  z 


(6b) 


F2^x'  y'  =  2  f  — 2~ — 7 — TTtl 

J0  +y  +z  LY'L 


z  x 

~2  1  ;  2x  2  2.1/2 

'  +  z  (x  +  y  +  z  )  ' 


(6c) 


This  solution  is  seen  to  satisfy  the  Laplace  equation  and  to  have  the  fol¬ 
lowing  properties: 

1.  Across  the  trailing  vortex  sheet  |  y - r| 1  <  6q/2, 
x  >  0,  z  =  0,  the  potential  6 has  a  jump 

6a>(x,y,  0+)  -  6<p(x,y,0-)  =  and  6<p(x,  y,0)  =  0 

elsewhere. 

2.  Across  the  trailing  vortex  sheet  6 cp^,  bcp ^  are 

continuous  and  bcp  may  be  discontinuous,  but 
1 

(6 (Py)  is  continuous. 

2  2  2 

3.  bcp  0  as  x  +  y  +  z  —  oo  except  over  the 
vortex  sheet. 

Here  bcp  is  conveniently  decomposed  into  two  parts,  bcp ^  and 
hep 2"  The  part  bcp ^  is  independent  of  x,  and  will  be  called  the  two-di¬ 
mensional  part  of  the  potential.  The  part  bep^  is  a  function  which  is 
odd  in  x  and  z,  but  even  in  (y-f]),  and  will  be  called  the  three-dimen¬ 
sional  part  of  the  potential.  The  whole  part  bcp  may  be  regarded  as  the 
fundamental  solution  of  a  vortex  element,  and  is  useful  in  constructing 
the  solution  of  more  complicated  boundary  value  problems.  By  using 
this  fundamental  solution,  the  velocity  field  generated  by  a  distribution 
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of  circulation  r(y)  along  the  wing  span  (-b^<  y  <  b^)  under  the  influence 
of  a  single  jet  can  readily  be  calculated;  this  is  given  in  the  next  section. 
This  result  will  be  further  generalized  to  the  case  of  multiple  jets. 
Finally,  the  approximate  lifting  surface  theory  of  Weissinger  will  be 
applied  by  taking  one  or  more  of  such  lifting  lines  appropriately  dis¬ 
tributed  over  the  wing  surface  to  represent  the  wing  together  with  ap¬ 
propriate  boundary  conditions  to  ensure  that  the  wing  remains  a  rigid 
surface. 

fl.  SINGLE  JET  THEORY 

1,  The  Velocity  Potential  of  a  Lifting  Line  Extending  through  a 

Single  Jet 

We  shall  first  derive  from  the  above  fundamental  solution  Eq.  (6) 
the  velocity  potential  of  an  infinitesimal  bound  vortex  of  strength  nqjSq 
placed  at  (0,  q,  0)  in  the  presence  of  a  single  jet  with  its  boundary  at 
r  =  1.  This  solution  may  further  be  regarded  as  the  fundamental  solution, 
or  the  Green's  function,  of  the  wing-jet  interference  problem.  It  turns 
out  that  the  determination  of  this  Green's  function  depends  on  whether 
the  vortex  element  is  located  inside  or  outside  the  jet  slipstream.  The 
velocity  potential  of  a  whole  lifting  line  extending  through  the  jet  can 
then  be  obtained  by  superposition  of  this  Green's  function.  We  shall  con¬ 
sider  the  two-dimensional  and  the  three-dimensional  part  of  the  solution 
separately. 


a.  Two-dimensional  Part  of  the  Velocity  Potential 
The  two-dimensional  part  of  the  velocity  potential  can  be 
readily  derived  from  Eq,  (6)  by  making  use  of  the  simple  reflection  into 
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(5) 

the  circular  jet  boundary,  as  given  in  Koning's  theory'  .  The  following 
two  different  cases  must  be  distinguished. 

(1)  1 1|  |  >1,  the  Vortex  Element  is  located  Outside 
the  Slipstream 

In  this  case  it  can  be  shown  that  the  velocity  potential 

outside  the  jet  is  due  to  the  original  vortex  of  strength  r(Tl)6q  at  y  =  r| 

2  2-1-2 

plus  an  image  vortex  of  strength  (1  -  p  )(1  +  p  )  r)  r("n) 6t|  at  y  =  1  /r\ 
(which  is  inside  the  jet).  On  the  other  hand,  the  velocity  potential  inside 
the  jet  is  due  to  an  equivalent  vortex  of  strength  2p(l  +  p^)  ^r("n)6q  located 
at  the  original  position  y  =  rp  For  brevity  we  write 

q  =U  -  H-2)/<l  +  lA  t2  =<1  -  |i)2/(l  +  n2).  (7> 

which  also  satisfy  the  relations 

(1  -  e2)  =  (1  -  tj)/p  =  p(l  +  £j).  (7a) 

Then 

1  ^  i  1 

b(pl  =  -p:  r(ri)6ri{F1(y-71,  z)  4  Fj(y-  z)}  for  r  >  1, 

(8) 

= -^r(T])6ri(l  -  £2)F1(y-n,  z)  for  r  <  1, 

where  r  =  (y^+  z^)^^  and  F^(y,  z)  is  given  by  (6b).  That  this  solution 
satisfies  conditions  (4)  and  (5)  is  readily  verified  by  direct  substitution. 

In  the  above  expression  the  terms  with  and  6^  represent  the  effect 
of  the  slipstream  on  5<^  outside  and  inside  the  jet  respectively.  If  the 
value  of  Th)  is  given,  then  it  is  seen  at  once  from  (8)  that  for  a  vortex 
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element  located  outside  the  jet,  the  effect  of  the  slipstream  is  to  increase 
<pj  outside  the  jet  and  to  decrease  (p^  inside  the  jet. 


(2)  |  r{  i  <1,  the  Vortex  Element  is  located  Inside  the 

Slipstream 

In  this  case  it  can  also  be  shown  that  the  velocity 
potential  outside  the  jet  is  due  to  an  equivalent  vortex  of  strength  (1  -  e ^ ) T( "n) 6t| 
located  at  y  =  r),  and  the  potential  inside  the  jet  is  due  to  the  original  vortex 
r(r))6t)  at  y  =  t)  plus  an  image  vortex  of  strength  -6 jT]"^r(Ti)6r|  located 
at  y  =  l/t]  (which  is  outside  the  jet).  That  is 


6<?1  ~Tn  r(rl)6r)(1"e  2^F i(y"rl»  z) 

1  £  1  1 
=  ■5^r(Tj)6ri{K1(y-Ti,  z)-  -jF^y-  -,  z)} 

n 


for  r  >  1, 

(9) 

for  r  <  1. 


This  solution  also  satisfies  conditions  (4)  and  (5),  as  is  readily  verified. 
The  effect  of  the  jet  slipstream  is  now  to  decrease  the  value  of  cp ^  both 
outside  and  inside  the  jet  for  fixed  I\ 

The  system  of  equations  (8)  and  (9)  represents  the  two-dimension¬ 
al  part  of  the  Green's  function  of  the  wing-jet  problem.  For  a  lifting 
vortex  line  of  a  given  circulation  distribution  T(y)  extending  from 
y  =  -b^  to  y  =  b^,  the  two-dimensional  part  of  the  velocity  potential 
can  be  obtained,  as  the  problem  is  evidently  linear,  by  integrating  (8) 
and  (9)  along  the  lifting  line  with  the  contributions  in  different  regions 
kept  in  order.  Thus  we  have 


(10a) 
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HTiJFjty-Ti,  z)dr| 


(J  +  ^ntiJF^y-ri,  z)dri 


-  e 


for  r  <  1.  (10b) 


b.  Three-dimensional  Part  of  the  Velocity  Potential 
It  has  been  found  that  the  three-dimensional  part  of  the  velocity 
potential  in  the  presence  of  a  jet  slipstream  cannot  be  expressed  in  terms 
of  the  simple  image  system  alone  as  for  the  two-dimensional  part.  How¬ 
ever,  by  introducing  the  same  image  system  as  that  of  (p^  to  one  part 
of  cp ^  and  by  leaving  the  remaining  part  of  cp^  to  be  determined,  the 
analysis  can  be  somewhat  simplified.  We  proceed  as  follows, 


1.  r\ ,  >1,  with  the  Vortex  Element  located  Outside 


the  Jet 

In  this  case  we  write  in  the  form  which  contains 
partly  the  same  image  system  as  that  of  (8)  as  follows 

1  61  1  ' 
b(Pz  =T5r  rh)5r|{F2(x,  y-t],  z)  +  F2(x,  y-  z)+  cpj 


for  r  >  1,  (11a) 

=  ^.r(Ti)6ri{(l-e2)F2(x,  y-rj,  z)+  (p.  }  for  r  <  1,  (lib) 

r  r 

in  which  F2  is  defined  in  (6c)  and  pQ  and  cp j  satisfy  the  Laplace 

equation  and  are  to  be  determined  by  using  conditions  (4)  and  (5).  For 

this  purpose  it  is  convenient  to  write  the  Laplace  equation  in  polar 

Z  Z  l/ z 

coordinates  (r,8,  x)  with  r  =  (y  +  z  )'  and  0  =  arctan  (z/y),  so  that 


7 ¥*7-^ 

or  r 


2  1 
9  (p 

—  2 

90 


9JL  =0 

9x^ 


(12) 
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The  above  integrals  can  easily  be  integrated.  By  doing  so,  one  may 

i 

readily  see  that  conditions  (13)  and  (14)  show  that  cp  has  the  following 
behavior 

1  1  -  Z 

cp  =  0(  I  x  I )  as  x  —  0,  and  cp  -  0(x  )  as  I  x  I  -*•  oo. 

Hence  cp  and  its  derivatives  will  be  small  both  near  x  =  0  and  |x|  =  oo. 
This  result  exhibits  the  advantage  of  expressing  cp  ^  the  form  (11), 
with  the  simple  image  system  first  singled  out.  Furthermore,  (13)  and 

i 

(14)  indicate  that  cp  will  be  a  function  odd  in  both  x  and  0.  It  therefore 

i 

follows  that  cp  must  be  of  the  following  form,  which  is  readily  verified 
to  be  a  solution  of  (12), 

00  oc 

cp  =  cp Q  =  j  ^  sin  n0  y  S^(k,  r|)Kn(kr)  sinkx  dk  for  r  >  1,  (15a) 

n=l  0 

co 

=  cpj  =  ~  ^  sin  n0  ^  sj^  (k,  r))In(kr)  sin  kx  dk  for  r  <  1,  (15b) 

n=l  ^ 
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where  I  and  K  are  the  modified  Bessel  functions  of  the  first  and 
n  n 


s(o) 


:(j) 


second  kind,  S'  and  S'J/  are  arbitrary  functions  of  k  and  rj.  The 


i 

functions  K  are  chosen  for  r  >  1  and  I  for  r  <  1  so  that  m 
n  n  Y  o 


and 


cp .  will  be  regular  respectively  at  r  =  oo  and  r  =  0.  The  super-indices 

J 

(o)  and  (j)  of  S  denote  the  regions  (outside  or  inside  the  jet)  with  which 
S  is  associated;  this  notation  is  convenient  for  extending  the  analysis  to 
the  case  of  multiple  jets. 

To  facilitate  the  determination  of  and  it  is  convenient 

n  n 

to  express  conditions  (13)  and  (14)  also  in  terms  of  Fourier-Bessel  inte¬ 
grals.  By  using  the  relations  (which  can  be  deduced  from  the  results  on 
p.  388  of  Ref.  15) 


(x2+T^-2r|  cos  0  +1)'^ 


oo 

■y. 


k(t)  -2t)  cos  8 


+  1)1//21 


sin  (kx)kdk, 
(16a) 


T)  IX 


x  (£2+t]2-2t]  cos  0  +  l)^2 


I  f  Sm.  kx  dk  f  tK  [  t(r)2-2r|  cos  0  4  1)1//2]  dt,  (16b) 

^  Jo  k  Jk/ 1  -n  J  0 


and  the  addition  theorem  (see  Ref.  15,  p.  361) 


sin  0  Kq[  t(rj2-2t]  cos  0  +  l)^2] 
oo 

=  ^  [  In_i(t)Kn_i(Tlt)"  In+1^t^Kn+l^t^  sinn0  for  |  r)  I  >  1,  (17a) 
n=l 


co 

=  ^  [  W^n-l^-  In+l<Tlt,Kn  +l(t)]  sin  nG  for  1  ^ 1  <  (17b) 

n=l 
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together  with  the  integral  (which  can  be  verified  by  differentiation) 


2n 


IE  n 


IJat)Kn(bt), 


(18) 


one  finds  by  straightforward  calculation  that  (13)  and  (14)  can  be  written 


(<p 


j'Wo> 

J  r  =1 


4  t  00  .oo 

-4rl |19) 


n=l 


i  i 


'  9<p.  9^o 

"ET, 


r=l 


00 

46 ,  v-  r»oo  ,  ,  , 

—  ^nsinnGj  [  y  -  )Kn(k)  +  In(k)Kn(kr|)]  sin  kx  dk,  (20) 


n=l 


i  i 

wnere  I  and  K  denote  the  derivatives  of  I  and  K  with  respect  to 
n  n  n  n  r 

their  arguments, 

A  remark  should  be  made  here  with  regard  to  the  argument  of  the 
Bessel  functions  when  it  becomes  negative.  Since  the  right-hand  sides 
of  (13),  (14)  and  the  left  side  of  (17)  are  not  changed  if  r|  is  replaced 
by  (-r|)  and  0  by  (tt-9),  the  right  side  of  (17a)  and  (17b)  is  also  unchanged 
under  the  same  transformation.  Therefore  we  may  adopt  for  the  purpose 
of  later  calculation  the  conventions 


Kn(-z)  =  (-)nKn(z),  In(-z)  =  (-)nIn(z),  (21) 

which  are  consistent  at  least  for  the  present  analysis. 
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Now  application  of  conditions  (19)  and  (20)  to  (15)  yields 

la\  y  f  2kl  i'k  (k-n)  ,  > 

1,  -  !  (— , \  , 

1  h-e  k(l  K  +K  i  )  "  i  J 

Inn  n  n' 

m  [i  K  +K  I*  ]  K  (krj) 

S^>(k,  ,,)  =  £  (1-6  )h  -JL.il  V"  ",  . 

1  2  11  l-e,k(i  K  +  K  I  ) 

1  '  n  n  n  n' 


(22a) 


(22b) 


In  these  equations,  for  those  Bessel  functions  whose  argument  is  k, 
the  argument  is  omitted  for  brevity.  Thus,  Eqs.  (11),  (15)  and  (22) 
constitute  the  three-dimensional  part  of  the  Green's  function  for  a  vortex 
element  located  outside  the  jet. 

(2)  |  t|  |  <1,  with  the  Vortex  Element  located  Inside 

the  Jet 

The  solution  of  this  case  can  be  obtained  by  the  same 
method  as  given  in  the  previous  section.  The  final  result,  as  can  be  veri¬ 
fied,  is 

6<?2  =  Tt?  r(Tl)6'rl {(l-fc  Z)F  z(xr  y_rl’  z)  +  ^o  ^  for  r  >  1,  (23a) 


6(p2  T^)61!  {F2(x,  y-ri,  z)  - F2(x,  y- -j-,  z)  +  }  for  r  <  1;  (23b) 


<P  o 


00 


2 
7 r 


sin  n0 


n=l 


oo 


n=l 


sin  n0 


T^(k,ri)K  (kr)sin  kx  dk 


T^(k,  ri)In(kr)  sin  kx  dk 


for  r  >  1  (24a) 


for  r  <  1;  (24b) 


and 
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(25a) 


tJVi)  =  6l<1-£2) 


,  [I  K  4K  I  ]l  (krj) 
2n  n  n  n  n  n'  1 


l-£,k(l  K  +K  I  ) 
1  '  n  n  n  n' 


m  2n  f  2kK„K„ijkn) 


n  n  n 
T 


-e,k(l  K  +  K  I  ) 

1  '  n  n  n  n 


K  (-))  . 

n'  r|  j 


(25b) 


In  (25)  those  Bessel  functions  whose  arguments  are  omitted  have  the 
argument  k. 

For  a  whole  vortex  line  of  distribution  T(y)  ranging  from  y=  -b^ 
to  b^  and  extending  through  a  jet,  the  solution  of  <p  ^  is  obtained  im¬ 
mediately  by  integrating  the  Green's  functions  of  part  (1)  and  (2).  This 
expression,  however,  will  not  be  explicitly  written  down  here,  but  will 
be  included  in  the  final  form  of  the  total  perturbation  potential  given  as 
follows.  In  summary, 


<p(x,  y,  z)  =  ^(y,  z)  +  cpz(x,  y,  z) 


(26a) 


and  with  the  notation  <p  =  <p  for  r  >  1  and  cp  =  cp,  for  r  <  1,  then 

J 

4 Vcpo=  f  ^  T(r])F(x,  y-Ti,  z)drj  -  £2j  ^  F(r|)F(x,  y-r|,  z)dr| 


+£l(L/^1'?F(x’y^■z)d,' 


00 


sinn8\  K  (kr)sinkxdk 
)0  n 


^(lOsJ^k,  T))dr| 


+  J  Ti)dri  , 


(26b) 
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.  J 


b 

-b 


2  f(Ti)F{x,  y-ii,  z)dri  -  £  2  ^  J  +  J  2  J  r(ri)F(x,  y-rj,  z)dt] 


1  00  oo 

-  1^- F{x,  y- I  ,  z)dTi  +  sin  n0  J  In(kr)sinkx dk  J 

“  1  T1  -  0  •  “ 


n=l 


r  n-1 
b 


r 


+  \  2  r(T1)s[ij)(k,  t, )dr,  +  j*1  r(T!)T{nj)(k,  TiJdT, 


(26c) 


where 


F(x,  y,  z)  =  Fj(y,  z)  +  F2(x,  y,  z)  =  ^ 

y  +z 


1  4 


:  2  2  2,1/2 

(x  +y  \z  ) ' 


.  (26d) 


By  using  the  rule  (21),  it  may  be  noted  from  (22)  and  (25)  that 


s[i°'j)(kI-T1)  =(-)n+1s[i°,j)(k,ri),  T^0*  j)(k,-n)  =(-)n+1T|i°’j)(k,T1) 


Hence,  if  b^  =  and  if  T(-r\)  -  T(rj ),  then  the  even  terms  in  the  infinite 
series  all  cancel  out;  similarly,  the  odd  terms  vanish  if  Tf-tl)  =  - r("n ). 

It  therefore  follows  that,  with  =  b^, 


<p(x,  -y,  z)  =  j-  (p(x,  y,  z)  for  f( -rj )  =  +  r(r|). 


As  Vj -*  Vq,  or  [j,  -*■  1,  we  have  £j,  £.>  ®  (see  Eq.  ?)»  then  only 

the  first  term  on  the  right  side  of  (26b,  c)  remains,  which  is  known  to  be 

the  lifting  line  potential  in  an  otherwise  uniform  flow.  On  the  other  hand, 

if  Vq  =  0,  V.  >  0,  then  £^  =  6^  =  1;  thus  (26)  reduces  to  the  potential  of 
J 

a  lifting  line  extending  through  a  free  jet,  that  is,  for  r  <  1, 
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The  solution  for  r  >  1  of  course  lose6  its  meaning.  In  the  general 

case  the  terms  with  the  factors  and  t  represent  the  effect  of  the 
jet  on  the  potential  <p. 


2.  The  Fourier  Expansion  for  Circulation;  Lift  Distribution 
The  spanwise  lift  distribution,  f( y),  defined  to  be  the  local  lift 
on  a  unit  span  of  the  wing  at  station  y,  is  given  by  the  Joukowsky  theorem 


At  the  boundary  of  the  jet  (y  =  4  1),  the  lift  distribution  must  be  continuous. 
Hence  it  follows  from  (28)  that  the  conditions 

r(l-O)  =  pr(l+0)  and  n-1+0)  =  ixr(-l-O)  (30) 

must  be  satisfied.  These  conditions  show  that  the  circulation  T  will 
have  a  jump  at  the  jet  boundary.  It  should  be  pointed  out  here  that  the 
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continuity  of  the  wing  surface  at  the  jet  boundary  requires  further  that 
the  downwash  be  continuous  at  the  boundary;  this  condition  remains  to 
be  enforced,  (see  Equations  45,  46). 

In  the  following  sections  dealing  with  the  single  jet  problem,  we 
shall  limit  ourselves  to  the  case  bj  =  b^  =  b,  i.  e. ,  a  single  jet  centered 
at  the  mid-span.  It  is  convenient  to  divide  the  span  into  two  regions: 

R  :  1  <  ly!  <  b  and  R.:  I  y  I  <  1, 

o'  1  7 

and  define  the  angles  iJjq,  4^  (see  Fig.  2)  by 

y  =  b  cos  ijj  for  1  <  I  y  I  <  b,  ( 31a) 

y  =  b  cos  =  cos  for  ly!  <1,  (31b) 

so  that 

4>o  =  cos  ^  =  p  and  ipj  =  0  at  y  =  1.  (31c) 

In  order  to  account  for  the  possible  slipstream  rotation  of  small 
magnitude,  the  circulation  T(y)  is  not  to  be  limited  to  even  functions  of 
y.  We  now  assume  that  T(y)  can  be  represented  by  the  following  Fourier 
series: 


T(y)  =  4V.rv(y)  in  region  R v>  v  =  0,  1,  (32a) 

2N 

o 

-  o(y)  =  ^  sin  mPo  in  Rq,  (32b) 

n=l 

2Nj 

ri(y)  =  ro(y)  +  ^  sin  mJjj  in  Rj  (32c) 

n=l 
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in  which  the  numbers  Nq  and  Nj  may  be  set  as  large  as  we  wish,  or 
as  practical.  The  expansion  for  T Q( y )  is  suggested  by  the  usual  lifting 
line  theory.  In  (32c)  the  constant  is  introduced  in  order  to  satisfy 

the  jump  conditions  (30),  the  difference  between  T^(y)  and  T0( y )  +  ^ 

in  lyl  <  1  is  chosen  to  be  a  sine  series  since  first,  the  half-period  Fourier 
expansion  in  0  <  iJjj  <  ir  is  rather  arbitrary  and  second,  the  sine  series 
is  more  amenable  to  analysis  than  the  cosine  expansion.  Now  application 
of  condition  (30)  to  (32)  yields 


N  -1 
o 

A(q1}  =  (p  -  1)^>  A^sin  (2n+l)p,  (33a) 

n=0 
N 

o 

0=^  4°J  sin  2nf3  .  (33b) 

n=0 

If  there  is  no  slipstream  rotation  and  if  the  angle  of  attack  is  symmetric 
about  the  mid-span,  then  A^  =  A ^  =  0  and  hence  (33b)  is  automatically 
fulfilled. 

Substituting  (32)  into  (29)  and  integrating,  we  obtain 


L 


8pVf 


N  -1 
o 


+  (fx-l)b 


k(°) 

2n-fl 


n=0 


sin  2n(3  sin  2(n+l)|3 
4n  '  — 4{n+T)  ' 


(34) 


It  is  noted  here  that  the  terms  with  A^  for  n>  1  and  aL^  make  no 

n  2n 

contribution  to  the  total  lift. 
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3 ,  Downwash  Distribution 

(8) 

According  to  the  approximate  lifting  surface  theory  of  Weissinger  , 
the  wing  and  its  wake  are  represented  by  a  bound  vortex  line  located  at 
the  quarter  chord  (which  is  also  the  center  of  pressure)  together  with  a 
system  of  trailing  vortices  as  in  the  lifting  line  theory,  and  the  flow  is 
then  required  to  be  tangential  to  the  wing  surface  at  the  three-quarter 
chord  points.  Let  the  quarter  chord  be  located  at  x  =  0  and  let  c(y) 
be  the  half-chord  length  at  y,  then  the  downwash  at  the  three-quarter 
chord  is 

w(y)  =  -8 cp(c,  y,  0)/8z  .  (35) 

In  carrying  out  the  differentiation  we  note  from  (26d)  that 

fz(C,  °)  =_L^{U — j-rryT.}  <36a> 

(y-n)  L  Uy-n)  +c  J  ■> 

"7  Fz(c’  y"  V  0)  =  ^V^]G(y'  V  c)  =  "  ^  G(y' c)  ’  (36b) 

•q  q 

where 

G(y,  c)  =  G^y)  +  G2(y,  c),  (37a) 

,2  1/2  2  1/2 

Gl(y)=f  G2(y,c)  =  ±[  (1+^)  -l]=-^/[l+(l+^)  ].  (37b) 

y  y  c  c  c 

The  function  G(y,  c)  is  decomposed  as  shown  above  such  that  G^fy,  c) 
is  regular  for  all  y  and  for  c  >  0,  the  only  singular  part  of  G  being 
G/y).  It  can  be  seen  that  G/y)  actually  corresponds  to  the  value  of  G 
at  large  distances  behind  the  wing  (the  so-called  Treffetz  plane  solution). 
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By  using  (36),  the  integrals  containing  F  in  the  expression  for  w  can 
be  integrated  once  by  parts.  Doing  so,  we  find  that  the  terms  evaluated 
at  the  limits  y  =  +  1  and  y  =  -  b,  b  all  vanish  on  account  of  condition 
(30),  relations  (7a),  and  the  condition  that  p(-  b)  =  H  b  )  =  0.  Furthermore, 

H  l  Kn(kr).in  ne]  z=0  =  1  Kn(kr)ain  „e]  z=0  =  f  Kn(ky) 


which  is  valid  for  both  y  >  0  (6  =  0)  and  y  <  0  (0  =  n)  provided  the 

continuation  of  Kn(kr)  from  0  =  0  to  0  =  7r  follows  the  rule  (21).  The 

same  result  applies  if  is  replaced  by  1^.  In  summary,  we  obtain 

the  downwash  w(y)  =  wQ(y)  for  1  <  lyl  <  b  and  w(y)  =  w.(y)  for  J  y  I  <  1, 

J 

as  follows 


oo 


7r  i~j  y 

V  =1 


OO 


[f+f 

LJ-b 


K  (ky)sin  kc  dk 


S(v0)(k,  r))r(T])dr| 


i; 


T(v0)(k,T!)r(ri)  dri 


dr 


(38a) 


4?TWj(y)  =  j  G(y-r|,  c)^  dr\  "  6  2 


I,! + 


+  e 


i£1G(y'Vc,^dT|  7  I, 


Iv(ky)sinkc  dk 


v=l 


rr 

•  -b 


^S(f{K  TiJntlJdT,  +  y  ^^(k,  T,)r( T|)dr, 


(38b) 


25 


Now  we  substitute  the  series  expansion  for  P  given  by  (32)  into 
(38),  we  note  that  some  of  the  resulting  integrals  can  be  evaluated  in 
closed  form.  For  example, 


p  cos  n0  d0  ^  sin  m[/ 

Jn  cos  0  -  cos  ip  ~  '  sin  i|j 


*  -1=6"  y)n+1”[  I yl  -(y2-l)1/2l  7(y2-l)1/2.  for  lyl  >  1, 

where  P  denotes  the  Cauchy  principal  value  of  the  integral,  and  (y^-lp^ 
represents  the  positive  branch  of  that  function.  After  some  rearrangement 
we  obtain  the  following  expression  for  the  downwash  w(y): 
a.  for  1  <  y  <  b 

2No 

=1  EAL0,{2^r+t2„(r'r;R  o+Ri>-£2tp„(ElRi,+  nnif  riRi>l 


e1b2[Qn(by;Ro)  +  An(by;^  ;RQ)] 


+(!-62)^  nA^  |^-2(sgny)n+1 


(y  -1) 


/  —  \  K  (ky)sinkcdk 

7  L  yJ0  v 

V=1 


b,  for  |y  <1 


{]  b+f  st0)(k,„ro(,)d,+yiT(°,(k,T1)r1(,)d,}; 


(39a) 


w.  r-'  /A\  f  sin  ndi 

>  =y  gA(o)-2 ° 

V7  /.on  sin  di 

J  „  - 1  L  0 


r +  n„(  r  ;Ro«i>- e  2I  p„<r  -Ro»+R„(^  -c  *„>! 


ib2[Q„(by;Ri)  +  An(byi^;R1)l} 


v  (i)  r  sinmK  i  r 

Z  ^l2^^>c;Ro+Ri)+£i 

n=l  L 


(1-y2)1/ 2  [  l+(l-y2r  2]  n 
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03 

7  5*0  lv(ky,si"kcdk{I. 

V=1 

fb  s(j)(k,  r| )  ro(n)dq  +  f  1  T(vj)  (k,  n^ri)  dr\  . 

Jj  J-1  J 


(39b) 


In  the  above  equations  the  new  functions  P^,  Q^,  f2n  and  A^  are  defined 


as  follows: 


2  pP  f* 


0 

">7r  ^  cos  n0  d0 


cos  n0  d0 
cos  0  -  £,  ' 


Pn^’Rl^  7T  J  cos  0  -  t, 


(40a) 


(40b) 


2  rp  r ' 


cos  0  cos  n0 
1  -  £  cos  0 


_  ,  „  .  2  (l7T"P  cos  0  cos  n0  /aik\ 

Qnlt;Rl>=  »Jp  1  -  rcos-IT  d9  ’  l4IM 

r(3  nt/Z 

Si  U;X;R)  =  f  U;X;0)d0,  R  U;X;R,)  =  fU;X;e)d0,  (42a) 

U  °  «Jq  n  11  1  Jp  11 

t  <co=  4(-)ntV°B6n)<:°;^ 'I  (42b) 

n  U+[l+X2(cos  0-  r)2] 1/2  1+  [1+ \2(cos  0  +  r)2]  1/2J 


n  "  (l+[in"(cos0-UlV6  1 

fy^XjR^)  =  S2nU;X;  R0)+  R^^jXjRj). 


An(^;Ro)=J  gn(C^;6)de,  An(&;MR1)=J  gn(£;X;0)d0,  (43a) 

o  p 

/r.x.oi-  k2/  (sec  6  -  £)cos  n6  ,  .n+1  (sec  0  +  r  )cos  n8  1 

8"a'M,--ii+liu2(sece-u¥/2  ('’  i+[ h\2(s.c e  +  y2] 1,  2  J ' 


(43b) 


An(^;\;R0+Rl)  =  Ah(^;X;R0)  +  AJ^XjRj). 


The  functions  Rn(  X)  and  A^U*  X)  are  readily  seen  to  be  regular  functions 
for  finite  f,  and  X  (that  is,  for  c  >  0).  The  functions  and  Q  defined 
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by  the  above  integral  representation  can  be  integrated  in  closed  form, 
and  P^,  (for  n  >  1)  can  be  expressed  in  terms  of  elementary  functions 
and  Pj  and  Q^.  For  the  purpose  of  numerical  computation,  however,  it 
is  quite  convenient  to  apply  the  modified  Simpson's  rule  directly  to  evalu¬ 
ate  these  integrals.  For  large  values  of  n,  the  most  significant  contri¬ 
bution  from  the  integrand  comes  from  the  neighborhood  of  the  points  of 
stationary  phase.  As  to  the  analytic  behavior  of  P^  and  Q^,  it  is 
noted  that  they  all  have,  a  logarithmic  similarity  at  the  jet  boundary 
y  =  +  1.  In  fact,  it  can  be  shown  that  if  £  is  a  small  positive  quantity, 
then  for  y  =  1  +  fc, 

Pn<fRl»£  t^HT1°8£  +  °M.  Qn(by;Ro)=54ilLl06£  +  0(1)  ; 

r  7Tb  sin  p 

(44a) 

and  for  y  =  1  -  t , 

P„<E;Ro'a- Wl0ge  +  °(1)-  +  °<1>- 

r  7Tb  6in  p 

(44b) 

The  singular  behavior  of  these  functions  at  y  =  -1  can  be  deduced  from 
the  above  results  by  using  the  relations 

PnK;R)  =  Hn+1PnU;R) ,  QnK;R)  =  Hn+1Qn(£;R) .  (44c) 


Aside  from  the  logarithmic  singularities  of  P^  and  Q  ,  we 

further  note  from  (39)  that  the  downwash  w(y)  also  has  square  root 

2  -1/2 

singularities  at  the  jet  boundary,  which  are  proportional  to  (y  -1)  ' 

or  (1-y^)  as  they  appear  in  the  series  ^  nA^.  The  strengths  of 
these  square  root  singularities  are  seen  to  be  independent  of  both  the 
span  b  and  chord  c,  whereas  the  strengths  of  the  logarithmic  singu¬ 
larities  (see  Eq.  44)  depend  on  b,  but  not  on  c.  This  implies  that  all 
these  singularities  are  solely  due  to  the  effect  of  the  slipstream  boundary, 
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and  therefore  are  inherent  in  both  the  lifting  line  and  lifting  surface  theory. 

Now  we  must  ensure  that  the  downwash  does  not  become  infinite 
at  the  slipstream  boundary  if  the  wing  surface  is  continuous  there.  In 
order  that  the  square  root  singularities  of  w(y)  at  y  =  1  and  y  =  -1  shall 
be  eliminated,  the  relations 


I  ■ 0 


(45a) 


2n =  0 
2n 


(45b) 


must  therefore  be  satisfied.  Furthermore,  in  order  to  remove  the 
logarithmic  singularities  of  w(y)  at  y  =  +  1,  the  additional  relations 


N  -1 
o 


(2n+l)A^+1  cos  (2n+l)p  =  0 


(46a) 


2nA^r|  cos  2n[3  =  0 


(46b) 


must  also  be  satisfied.  Conditions  (45a,  b)  and  (46a,  b)  may  be  regarded 

as  four  constraints  on  the  coefficients  and  A^  so  that  the  down- 

n  n 

wash  will  be  finite  at  the  jet  boundary. 

It  is  convenient  for  the  purpose  of  numerical  computation  to  first 

remove  the  singularities  in  the  expressions  for  w  and  w.  by  applying 

o  J 

the  conditions  (45),  (46),  or,  which  is  equivalent,  by  making  in  (39a,  b) 
the  following  replacements: 
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-2(6gn  y) 


n+1  [ly|-(y2-I)l/2]n  „_|n+l  l-[  I  y  | -(y2-!)1/2]  n 


(y7^7 


by  2(sgn  y) 


(y7^ 


(i^ luuCw-  by  ' (8en y)"'1}’ 


P  (  X  ;R  )  by  P*(£  ;R  )  =-  f 
nv  b  o'  7  n'b  o'  tt  J 


p  r  1  +  (-1)"*1 


0  Lcos  0  - ^  cos  9  +  ^  J 


p„(e'ri>  by  p!(c'ri)  = 2 


n'E,rV  ~1 


,/2r  i 


(-i) 


n+1  -i 


cos  0  -  ^  cos  9+^ 


Qn(by;Ro)  by  Q>y;R0)  =| 


fp 

r 

cos  0 

,  >n+l  _  "] 

(-)  cos  0 

^0 

IT 

-by  cos  0 

+  1+by  cos  0 

Qn*by'Rl*  by  Qn<by;Rl*  =i  J 


7r/2 


a  /  »n+l 

cos  0  (-)  cos  0 

1  -  by  cos  0  T+by  cos  0 


(cos  n0  -  cos  n|3)d0  , 


(cos  n0  -  cos  n(3)d0, 


(cos  n0  -  cos  n|3)d0  , 


(cos  n0  -  cos  np)d0 


With  these  substitutions  the  downwash  w(y)  is  therefore  regular  every¬ 
where  over  the  wing.  Although  conditions  (45a,  b)  and  (46a,  b)  have 
been  once  applied  as  shown  in  the  above  substitutions,  they  must  still 

be  used  in  the  calculation  of  the  coefficients  and  A^  since  no 

n  n 

coefficient  is  actually  eliminated  in  this  process. 


4,  Application  of  Boundary  Conditions  on  the  Wing 
In  order  to  take  the  effect  of  small  slipstream  rotation  into  ac¬ 
count,  we  make  the  following  simplifying  assumptions, 

1,  The  axis  of  the  cylindrical  slipstream  is  parallel 
to  the  undisturbed  main  stream  and  passes  through 
the  mid-chord  point  of  the  wing.  The  velocity  of 
the  slipstream  may  have  a  rotational  component 
with  distribution  rco(r)  about  the  jet  axis,  where 
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Since  in  the  conditions  (45a,  b)  and  (46a,  b)  the  odd  and  even  terms  are 
not  coupled,  it  becomes  convenient  to  decompose  the  problem  into  even 
and  odd  parts  as  follows. 


31 


a.  The  Even  Part 


The  boundary  conditions  of  this  part  reads 

wj,  even(y)/Vj  =  a(y)  for  0  <  V  <  1  (49a) 

Wo,  even^y^Vo  =  for  1  <  Y  <  b  (49b) 

where  w.  and  w  are  respectively  the  even  part  of  w. 

and  w  ,  which  can  be  obtained  directly  from  (39)  by  deleting  the  even 
terms  with  n  =  2,  4,  6,  . .  .  ,  2N.  The  two  conditions  (49a,  b)  contain 
Nq+  Nj  unknown  coefficients  A^jj  ,  ,  which  should  be  deter¬ 

mined  under  the  constraints  (45a)  and  (46a).  Hence  we  may  take 
(Nq+  Nj-  2)  appropriate  points  in  0  <  y  <  b  at  which  we  apply  conditions 
(49a,  b).  For  instance,  we  may  choose  (N^-l)  points  in  0  <  y  <  1  and 
(Nq-1)  points  in  1  <  y  <  b. 

b.  The  Odd  Part 

For  this  part  we  have  the  boundary  conditions 

Wj,  odd^  =  y0J^  for  0  <  Y  <  1  ’  (50a) 

Wo,  odd(y^  =  °  for  1  <  y  <  b  ,  (50b) 

where  w^  and  wq  are  the  odd  part  of  Wj  and  wq,  which  are 

obtained  from  (39)  by  deleting  the  odd  terms  with  n  =  1,  3,  5,  .  .  .  .  These 
two  conditions  contain  (NQ+  N^j  unknown  coefficients  and  A^j  , 

which  must  also  satisfy  conditions  (45b)  and  (46b).  The  same  points 
chosen  for  the  even  part  may  also  be  used  for  this  part.  Thus  the  problem 
is  now  reduced  to  solving  a  system  of  (N  +  N^)  linear  algebraic  equations 
in  either  the  even  or  odd  part. 
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C.  MULTIPLE  JET  THEORY 

1.  Geometric  Configuration  of  a  Wing  and  Multiple  Jet  System 

We  suppose  that  there  are  in  general  J  pairs  of  jet  slipstreams 
which  are  symmetrically  located  with  respect  to  the  mid-span,  so  that  the 
total  number  of  jets  is  2J.  All  the  slipstreams  are  assumed  to  be  circular 
cylinders,  with  their  central  axes  parallel  to  the  main  stream  and  passing 
through  the  center-chord  points  of  the  wing.  The  span  of  the  wing  is  divided' 
into  (2J  +  1)  regions,  denoted  by  Ry,  v  =  0,  1,  2,  .  .  .  ,  2J,  which  are  defined 
as  follows  (see  Fig.  3). 

Rq  is  the  region  on  the  span  outboard  of  all  the  jets; 

R2v  j  is  the  region  inside  the  vth  jet  pair,  v  =  1,  2, .  .  .  ,  J, 
counted  from  the  wing  tip  toward  the  center; 

R^v  is  the  region  between  the  vth  and  the  ( v+l)th  jet  pair s. 

The  boundaries  of  these  regions  along  the  y-axis  are  marked  by  the  points 


■b0=  “b, -br -b2, . . -b2J,  b2J+1-  0,  b2J, . . . ,  b2,  bj,  bQ  =  b. 


For  every  point  y  in  the  region  R^  we  define  (v-tl)  angles,  denoted  by 


V  by 


y  =  bgC°s  S  =  0,  1,  2,  .  . .  ,  V,  for  bv+1  <  lyl  <  bv  .  (51) 

At  the  point  y  =  b^,  let 

ib  =  =  cos_1(b  /b  ),  s  =  0,  1,  2,  .  .  .  ,  v,  v=  0,  1,  .  . .  ,  2J  .  (52) 

s  s  vs 


Furthermore,  all  the  jets  are  assumed  to  have  the  same  radius,  which  is 
normalized  to  unity.  The  centers  of  the  vth  jet  pair  are  at  y  =  a^ 
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v  =  l,  2,  .  .  . ,  J. 


(53) 


and  y  =  -a^,  where 


1 

1 


(b2v-l+  b2v)f 


a  =  -a  . 
-  v  v* 


Finally,  the  geometric  incidence  of  the  wing  and  the  slipstream  rotation 
of  the  jet  pairs,  if  any,  are  both  assumed  to  be  small  and  symmetric  about 
the  mid-span. 


2,  Velocity  Potential  of  a  Wing  Extending  through  Multiple  Jets 
The  image  system  of  a  vortex  element  due  to  a  single  slipstream 
has  been  determined  in  Section  Bl;  it  will  be  called  the  primary  image 
system  for  the  problem  of  multiple  jet  slipstreams.  The  image  system 
of  a  vortex  element  generated  by  a  set  of  multiple  jets  contains  the  images 
of  the  primary  and  higher  orders  (i.  e.  the  images  of  the  primary  system 
and  in  turn  their  images)  since  it  takes  the  images  of  infinite  orders  to 
completely  satisfy  conditions  (4)  and  (5)  at  all  the  slipstream  boundaries. 
In  order  to  simplify  the  analysis,  the  image  system  of  a  vortex  element 
due  to  a  set  of  jets  will  be  determined  by  taking  into  account  only  the  pri¬ 
mary  and  secondary  images  and  by  neglecting  the  images  of  higher  orders. 
The  final  result  of  this  approximation  is  expected  to  be  accurate  enough 
for  practical  appUcations.  In  this  manner  the  two-dimensional  part  of  the 
velocity  potential  can  be  obtained  as  follows. 

a.  For  a  Vortex  Element  r("n) 6t|  located  at  y  =  T  in  R, 
(Outside  All  the  Jets) 


(1-ap 


cp  outside  all  the  jets  is  due  to  r("n)6T|  at  y  =  t|  plus 
1  1 

~2  r(r))6r|  at  y  =  a  + -  ,  summing  over  1,1  =  -  J,  ,  .  .  J; 


i 


cp .  .  inside  the  \th  jet  is  due  to  (l-e?)r(r|)6r|  at  y  =  r) 

‘,U-«2>  J'  1 

plus  - T-r(T))5r|  at  y=a  + - ,  summing  over  l ,  # 

Ol-a,)  ‘  ’1‘al 
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b.  For  a  Vortex  Element  r(~n) 6r|  located  at  y  =  t\  Inside 
the  \th  Jet 

( 0  outside  all  the  jets  is  due  to  (1-e ,)r(r|)6r|  at  y  =  r| 

€id-C2)  °  ^  2 

plus  - -j—  r( r|) 5 "p  at  y  =  a.  + -  ,  summing  over  1  ,  * 

<n-y  ‘  ^ i 


p.  .  in  the  <th  jet  is  due  to  (1- e7)*'T(r})br\  at  y  =  t| 
J  #  * 


for  i,  ^  X.; 


(p  .  .  in  the  X.th  jet  is  due  to  r(T])6r]  at  y  =  r), 

el  J'  1  el(1_c2>2  1 

plus - ^  r(T|)6r|  at  y  =  a  +  — -  ,  plus  - 2  r(rj)6ri  at  y  =  a 

(T|"a\.)  ^  ^'ai)  1 

summing  over  l,  $  X. 


Again,  like  the  previous  case  of  a  single  jet,  the  three-dimensional 
part  of  the  potential  can  be  decomposed  into  two  components,  of  v/hich  one 
part  can  be  written  down  according  to  the  same  image  system  as  given 
above,  and  the  rest  can  be  expanded  in  a  Fourier-Bessel  series,  Since 
it  has  been  found  that  the  contribution  of  the  Fourier-Bessel  series  is 


very  small,  this  part  of  the  velocity  potential  will  be  constructed  by 
including  only  the  images  in  the  immediate  neighboring  regions,  whereas 
the  images  in  the  other  regions  and  the  higher  order  images  will  be  neglected. 
Thus,  summing  up  all  the  vortex  elements  along  the  span,  we  obtain  the 


following  result  for  cp.  With  (p  -  cp  outside  all  the  jets  and  (p  =  cp.  . 

o  J,  a. 

inside  the  \th  jet  (which  will  be  taken  to  be  one  on  the  positive  y-axis), 


we  have 
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J  -b,  ,  b, 

nZvil  n  Zv  f  ^ 

+J  {f(x.  y-n.  *)*£ 


'•'-b  ^b  v 
v=0  °Zv  °Zvil 


t>y'ai-n^'z|} 


)|  r(tl)dr) 


J  “  b  -%  t)  » 

r~'  n  Zv  o  2v  - 1 


♦<M2,yr  v+f  v  {F(x,y-^Z)+y  — L^F(x,y-a  --i-tz)i 

v^l-b2v-l  Jb2v  1  £,l^S  1  J 


w  ‘Mi> 


r(ti)dri 


J  00 


b2i  rb2i-2 


|  2  ^  Sinm0ij1  Km(kri)sinkxdk  {J  +  J  Sm<k»  'l“ai)r(ri)dT1 


/  =-J  m=l 


b2i+l  b2i-l 


Tm  (k'T1"ai)r(Tl)dT1  ’ 


(54a) 


J  -b,  tl  b, 
r —  p  2v-fl  p  Zv 


L-J  Ou  * '  K 

v=0  °2v  °Zvil 


I  {F(x-y-’<’ 


f  ’  £1  !  1  r  2^-1 

)  - ,  F(x,  y-a  -  — — ,  z^I^dr]  +  \  F(x,  y-7i,  z)r(Tl)clri 

V.  (n-a  T  i  11_a;  J  Jb„ 


i=-j  (,1-ai 

(#M 


J  ,  b 


+  dV 


.  2v-l 


J  -b. 


F(x,  y-T),  z)fdT]  + 


F(x,  y-T],  z)rdT] 


Mr]  J 


v=l  2v  v=l  -b,  , 

(*M  2,M 

,b2\-ir  s  f  ei(1-ez)2  i 

(#M 


^a#  J 


)rr(T))dr] 


b2\  r.b2\-r]  ... 


+  4^  slnm0x  C  Im(krx)sinkxdk{£  +£  ’1sW(k,^-ax)rdn 


b2\+l  hZ\-l 


P  2X.-i 

U  Tm  (k^-aX)rdT] 
b2\ 


(54b) 
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where  in  the  second  term  of  (34a),  the  sign  )  denotes  the  sum  over  i 
except  when  £  =  v  for  p  ^  0  and  f  =  -v  for  q  <  0  (i.  e.  ,  the  v  th  jet 
pair  are  deleted  in  the  £ -summation.  In  the  above  equations  F(x,  y,  z) 
is  given  in  (26d)  and  S(o),  S(j),  T(o),  by  (22)  and  (25).  The 


m  rn  m  m 


variables  r  ,  8  are  defined  by 
f  i 


y  -  af  =  rfcos  V  z  =  rfcos  V  * 


(54c) 


3.  The  Fourier  Expansion  for  Circulation;  Lift  Distribution 


The  spanwise  lift  distribution  £(y)  and  the  total  lift  L  are  again 
given  by  (28)  and  (29).  The  circulation  distribution  F(y),  and  hence  the 
local  lift  f(y),  is  assumed  here  to  be  an  even  function  of  y,  that  is, 
r(-y)  =  T{ y )•  In  order  that  the  pressure  and  the  local  lift  shall  be  con¬ 
tinuous  at  the  jet  boundaries,  y  =  b.,  b^,  .  .  .  ,  b^j,  the  following  conditions 
(similar  to  30) 

f  P  0  r  v  =  1.  3,  .  .  .  2J-1, 

T(b-0)=\  ^Hb+0)  ior\  (55) 

V  ll/pJ  V  l  v  =  2,  4, .  .  .  2J  , 


must  therefore  be  satisfied.  The  continuity  of  the  downwash  at  the  jet 
boundaries  will  be  examined  later  in  the  analysis. 

The  circulation  T(y)  of  the  present  case  is  assumed  to  have  the 
following  Fourier  expansion 

r(y)  =  4V.rv(y)  ,  for  bv+1<lyl<bv,  v  =  0,  1,  2,  .  .  . ,  2J  ,  (56a) 

N  -1 
o 

ro<y)=^  A2n  +  i  sin  (2n+1H0’  (56b) 

n=0 
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N  -1 

v 


rv(y)  -  r^y)  +  A(oV,  +  \  A^jSin(2n+l)4ivl  v=l,  2. . 2J.  (56c) 


u-0 


The  last  two  equations  may  also  be  combined  to  give 

N  -1 

v  .  s 


r„(y)  =  ),  {  AoS)  +  ),  A2n+1  si»  (20+1)4,1  for  v=0,  1,  2 . 2J,  (56d) 

s=0  n-0 

=  0  being  understood.  The  above  expansion  gives  r(y)  an  even 
function  of  y,  as  clearly  suggested  by  the  fact  that  the  downwash  w(y) 
will  be  symmetrical  (even  in  the  presence  of  slipstream  rotation)  about 
the  plane  y  =  0. 

It  should  be  mentioned  here  that  a  few  different  expansions  for 

r(y)  have  been  examined  by  the  authors  by  performing  several  numerical 

programs  on  an  IBM-709  machine.  One  of  these  expansions  is  particularly 

worth  noting.  This  expansion  assumes  the  same  form  as  (56)  except  that 

in  (56c)  the  term  ry  ^(y)  is  removed.  This  reduces  the  number  of  terms 

of  the  expansion  in  the  inner  regions  R  ,  v?  1.  Thus  the  series  with 

( v ) 

coefficients  is  assigned  only  to  the  region  R^,  in  which 

bv+l  <  |7|  55  bv’  or  0  ^  +  and  ir~  P[,V^  <  Strictly 

speaking,  these  series  in  this  type  of  expansion  are  not  Fourier  series 
since  the  set  of  functions  sin  (2n+l)4jv  are  not  mutually  orthogonal  over 
the  region  R^.  It  has  been  found  that  this  type  of  expansion  induces 
numerical  instability  in  the  computing  program.  The  expansion  (56),  how¬ 
ever,  has  been  found  satisfactory. 

Now  application  of  conditions  (55)  to  (56)  yields 
v — 1 


A 


‘H(l  1  A»S,+i  AZnil  sin(2n+1,P<8V)l} 


(57) 


s=0  n=0 

v  =  1,  3,  .  .  . ,  2J-1, 


for 


{ 

tv  = 


2|  » p  2J  p 
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in  which  Aq  =  0.  Thus  the  total  number  of  the  arbitrary  coefficients 


^  *■ 


Ns  =  Nq  +  Ni  +  .  . .  +  N2j  . 


The  total  lift  L  is  obtained  by  integrating  i(y)  along  the  span. 


8p  V  4Vf  J0 

J  J 


Oq  D 

7T  j  Vlocal  r(y)dy  =  f1  I  £  rv(y)dy 


v  =  0,  1.  ..  v+1 


J  b. 


+  (1  -  H) 


r2v-l(y)dy 


v=l  2v 


On  substituting  (56d)  in  the  above  integrals,  and  making  use  of  b^  i"^2v=  2 


(which  is  the  diameter  of  the  jet),  we  obtain 


bviAiv,+4’Ai‘')]+<i-t'> 


J  Zv-l 


2A(s,+  *  A<s>[  p<2v>-  p(2v-1>] 
o  2  1  L  J 


V=1 s=0 


N  -1 

bsfV  (s)  (s)  sin  Znp[ZV}>  ~  sin  2npg2V_1* 

(T  Z  <Alni,-A^l)  - - - ZH - — 


(s)  sin  2N  p-2v^-  sin  2N  (3^2v_1^ 
srs  s  rs 

L2N  -1  2N 

s  s 
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4.  Downwash  Distribution 


Again  we  let  c(y)  be  the  half-chord  length  of  the  wing  at  station 
y.  Then  the  downwash  at  the  three-quarter  chord  can  be  derived  from 
the  solution  (54)  by  differentiation, 

w0(y)  =  y’  0)  f°r  hZUl<  lyU  b2\'  ^  = 

(60) 

0 

wj,  \M  1  •  Si  Q),  X(c' 01  for  b2\<  lyl®  b2\-r  b  =  1’2 . J' 

On  integration  by  parts  of  those  integrals  in  the  expression  for  d(p[c,  y,  0)/8z 
which  contain  F  ,  and  by  making  use  of  (36),  (37)  and  (7),  we  find  that 
those  terms  containing  G(y-q),when  evaluated  at  various  limits  of  integration, 
all  vanish  by  virtue  of  (55)  and  (7a),  and  we  finally  obtain 

ZJ  b 

( 1“ ^ v dti  e  Z )  j  [  G(y-ri,  c)  -  G(y  +  ri,  c)]  d-q 
v"=0  "bv  +  l 


(#v) 

+  4trw^  ,  (6la) 
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iz  l  l . 

4’Wo(',)  *  T  1  1  n»2v+0)[G(y.ai.E-ii-,c)  -  G(y+V  .  c)  ] 

(i*v) 


-6 


J  J 

2V  V 


i  Z  Z  r(b2v-i+°)  i  G<v  v  sr-ri-,  * c) '  G{y+a*  » c)  1 

v=l  t=- J  ^v_1  *  2v_1  * 


(i*v) 


2  V  V  f00  Km<k(y-ai>>  .  ,  „ 

¥  /  /.  m  -  sinkcdk 

i 


i=-J  m=l 


o  y-a. 


(I 


2/ 

2i+l 


I 


2i-2 


2/-1 


2/-1 


s|°^(k,  rj-a  )T  dt|  +  Z  T(o)(k,r, 

m'  i  ,,  m  '  1 


2.1 


-ai)rdT)|  ; 


2J  b 

47rwjf  x(y)  =  (1-£2)^  [G(y-n,  c)  -  G(y4ri,  c)]  -^Tdri 

v=0  bv+l 


(i*M 


+  4»Wj*x  ,  (61b) 
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v  1 

4'wj;x<y's-V1-ti>Z  ir(b2)i+0,G(',‘arB— 13-  -c> 


i=-J 

(i*M 


2\  i 


"fxr(b2\-l~f°)G(y-ai'  F^-~-a”  *C)] 


J  J  , 


+tl<1+tl>l  I  r(b2v40)[G(y-a(-  ‘  ,c)-G(y-ap  ' -  .0] 


V=1  f=-J 

(i*M  ) 
J  J  , 


2v  “i 


/  I  u  . 

2v  l 


>  2  r(b2v.140)[G(y-aJ-_i_a  , 

=1  f=-J  1 


i—  c)  -  G(y-a  +5-1-J-  c)] 

v-1  l  2v-\  i 


U*K  v) 


?  r00  UMy-aj)  (chzx 

)  m  \  -  sin  kc  dk  ■<  \ 


,b2k-2  rbZ\-l  ...  ■, 

Sm  <k.  1- aX>r  dl  +  \  (k,  r,-ax)  T  dn  L  . 

b2\-l  b2k  J 


In  (61a)  and  (61b),  the  notation  6°bb  is  defined  by  6°bb  =1  if  v  is  an 

* 

odd  integer  and  =0  if  v  is  an  even  integer.  In  the  expressions  for  wq 

*  -1 
and  w.  . ,  the  terms  with  G(y  +  a„  +  (b  -a.)  ,  c)  arise  from  the  fact  that 

the  reflections  of  the  image  systems  are  not  complete;  these  terms  will  be 
removed  if  all  higher  order  images  are  included  in  the  solution  of  (p.  How¬ 
ever,  it  is  noted  that  these  G  functions  are  regular  in  their  respective 
regions, and  their  magnitudes,  with  their  multiplication  factors  (depending 
on  p,  i  ].  or  6  accounted  for,  are  generally  very  small.  Furthermore, 
it  has  beer  found  from  the  numerical  results  that  the  contributions  of  the 
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)}t  >>  $ 

Fourier-Bessel  series  in  w  and  w,  .  are  al3o  small.  Hence  w 

o  j,  \  o 

$ 

and  w.  .  are  rather  unimportant  compared  with  the  remaining  terms 
J*  K 

in  the  expression  for  w  and  w.  . . 

°  Jr  X 

Now  from  (51)  and  { 5 6d)  we  find 

N  -1 
v  s 

=4Vj£  1  <2n+1>A2nilCO- (Zn+D+.d^  (62) 

8  g  =0  n=0 


for  bv+1<  ly  U  by, 

v  =  0,  1,  .  . . ,  2J. 


Using  this  expression  for  T  dr)  and  (36)  for  G,  we  find  that  some  of  the 
resulting  integrals  which  contain  G^(y-r))  and  G^(y+r])  can  be  evaluated 
in  closed  form.  The  final  result  is  given  in  the  following.  Since  w(y)  is 
evidently  an  even  function  of  y,  only  the  result  for  positive  y  will  be 
given  here. 


a.  Outside  the  Jets, 


bZ\+l< 


2X’ 


\  —  0plyt««y\J 


w0(y) 

v~ 

j 


=  2 


2N  -1  .  . 

2X  s  s 

nA  sin  mil 
n  Ts 


j  b  Sln 

s=0  n=l,  3, ...  fa 


2J 


-  2(1-6^) 

s=2X+l 


(!-6(s2ul)e2) 


2N  -1  .  , 

s  .(s) 


nA'  [  y  /b  - 

n  1  ’  '  b 


i 


(y/bJ2-i]n 


n=l,  3, 


(y/bs)  -  i 


T  ,  .  2N  -1  .  . 

J  2v-l  s  ^(s) 


i 

v=i  s=0  n=l,  3, 


-  (1-«S,U1,6i2U1>,[Pn,(f;R2v-1) 


s 


;R2v-1^ 
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‘"W”  (8) 

+  ‘2(i-6^)  ")  Z-2-  !  p(2U1) 


i^; 


v=2\+2 


t  p(»«)(.  y  .  y  Rw)1 


^  'V  '  *  •%  v  *\ 

v  =  2X.+ 2 

,,  2N  -1 

2J  vs  .  (s) 

*1  I 

v=0  s=0  n=l,  3,  ...  6  8 


+  «LS)  ( •  Rvn 


-‘i?  »-*?**  A 

v=0  f=-J 


2N  -1 
v  s 


v<dTlVZ)'I  7  -^U^vv 

s  =On=l,  3,  ,  .  . 


w!(y) 


+  Qj1S)(-y;a(;Rv)tAj13l(y;ai;RvMAjiS,(-y;a  ;Rv)]  +-^—  .  (63a) 


b.  Inside  the  Uh  Jet,  b^  <  y  <  b^  ^  \  =  I,  2,  .  . .  ,  J 

,  2N  -1  .  . 

cX.-l  s  (s 

v-  v-  nA  sin  mb 

=  2>  )  “ - A-8 


)  sin  ib 

s  Ts 


s=0  n=l,  3,  .  . . 


VS  [y/bs  -  s^/hs)2"1]11 


s=?A 


n=l,  3,  .  .  . 


(y/bs)  -  i 


2X-12Ns-1  (s)  Z\-Z 

>\5.  7 


n  'B~ 


s=0  n=l,  3,.  . 


■  t  Rv! 
v=2X. 


<2N ,,0-1)  (8) 


+  Pa8,<-f 'Z  Rv»-£2(1-£2)(1-6JX>  I 


V  a  rp(2U,y  .V 

■’f,  ^ir[  -  VA 

n=l,  3,...  v=2\+l 
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„  _  ,  2N  -1  .  . 

2J  J  2v  _s  (s) 


v=2\+l 


v  =  0  s  =0  n=l,  3,... 


*  [Pn8’(  f  •R2v)  +  Pn6’(-f  ;R2v)1 

2J  8  v  2Ns_1  8 

i  S 

v=0  s=0  n=l,  3,...  b  8 

2W  2Ng-l  ^ 

*  n‘6’(-f  i  T-  i  M  +  ‘f  1  1  nne’<f  •  T-  •  R2X-1> 

6  ,  ,  S  8 

s=0  n=l,3 

,,  ,  2N  -1 

Z\- 1  8  J 

-si  I  "A,.‘,,i{l,4Si<Z-e?»[4,’^l|>R2X-l» 

s=0  n=l,  3 

J  2y  2Ns_1 

+  A(11,,(T»iiR2x.1)]}-‘1(i-‘2)I  I  I  .“’I 

v=0  6=0  n=l,  3  i=-J 


*j.  & > 


Y.  • 


(63b) 


5?  v  m. 

In  the  above  equations,  wq  and  Wj  ^  are  given  in  (61);  6n  are  the 


Kronecker  deltas, 


6m  =  1  if  m  =  n  and  6m  =  0  if  m  £  n  ; 
n  n 


(63c) 


and  the  functions  P,  Q,  Q,,  A  are  defined  as 

o(v+i) 


p'*6’ (y/W  =  If 


cos  n0 d0 


6  -  y/b 


(64a) 


,(v+l) 


,  >  j  n^s  (b  cos  9  -  a  )cos  n0  d0 

Qis  (y;ai;Rv) '  4  \  i»i  - — — : 


rB  >Ds 

•Vv)  1  ♦ 


(bscos  0  -  ai)(a ^  y) 


(64b) 
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..  b  ,  b  .  n  (c  o  v/b  )cos  in  :9 

n(S>(>  :-£;R  |4(-S)  \  .  - - — - ,  (°-  I 

"  \  c  v  '  c  -yv)  u  2  2  1/2 

1  +[1  t  (  -|)  (cos  9  -  ^-)  ] 

8 


AnS)(y:ai:Rv)  =  ' 


P^V  +  J^(, - L - +  a  -  y)  cos  n9  d0 

vs  b  cob  0  -  a,  i  7 

s  i 


7TC  '  P 


1  +  (  1  +  -L  {, - - 

c  b  cof 

C  B 


t  ,2,  1/2 

+  a  -  y)  ] 


.  (64d) 


The  strongest  singularities  which  appear  in  the  expression  (63) 

for  w(y)  are  the  square  root  s ir.gularitie s  at  y  =  b  ,  s  =  1,  2,  .  .  .  ,  2J. 

s 

In  order  that  these  singularities  at  the  jet  boundaries  be  removed,  the 


following  2J  relations 


2N  -1 
s 


l 

n=l,  3,  .  .  . 


nA^°  =0,  s  =  1,  2,  ....  2J 


must  be  satisfied. 

•  ( 8 )  ( g ) 

Further,  it  is  noted  that  the  functions  P'  7  and  Q'  7  may  have 

n  n 

logarithmic  singularities  at  the  jet  boundary.  These  singularities  are 
listed  below. 

1.  In  y  ~  ^  —  1*  •  •  •  *  “l* 

(s)  .  b2X+l*fc  ,  )  —  b  0^S^b  +e-a  R  1 

n  (  5  ’  R2X+P  “  sUn  (b2k+l  +  €,a\+l’R2V 

9 


2  cos  n(3 

1  ~  TTk 

sin  pa 


2\+l 


log  e, 


(66a) 


Pn  ’  (-TT^  ;  R2X-1>  =  baQn ’<b2V  £;aX;R2X> 
s 


2  cos  npj 

f  ~  721 

sin  pa 


log  6  ; 


(66b) 
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2.  In  b?.  <y^  b?.  X  =  1,2,..,,J, 


PL8>I  -T“  !RZx'Sb.Q!,8><b2XWiVR2X-l» 

s 


,  cos  np‘ 

V  ~T2 !X 
sin  pc 


log  e  , 


(s)  .  b2\-l-£  .  ^.qIs)..  . 

n  (-T  ^ZX^  "  bsQn  (b2X-l  ,aX,R2X-r 
s 


(66c) 


^  cos  np* 
1 

sin  pc 


log  e  . 


(66d) 


To  remove  these  singularities  we  impose  the  following  conditions 


,,  ,  2N  -1  .  . 

2X-2  s  (s) 

c—  nA'  ' 


s=0  n=l, 3,... 


2X-1  2Ns-> 


n  04\-i  /  ,  nc.) 

g -  cos  nPs  /sin  p^ 

8 


0,  X  =  1,  2,  .  .  .  J  ;  (67a) 


^  cosnp^X/sinp^  =0,  X  =  1,  2,  .  .  .  J.  (67b) 

s=0  n=l,  3,.  .  .  6 


By  making  use  of  the  conditions  (65)  and  (67)  we  may  introduce  in  (63) 
the  following  substitutions: 


[£-  J(^)2-U7  (^)2-l  by  {[i-  J(*2-l  ]“-!)/, 


^7X4.1^  x  -  o,  i,  . . .  j , 


p(s)  i  y  \  uv  p^(®)  [  y  .  rj  \ 

n  'F-  ,R2X+1)  by  n  lF~  ’  R2X+1> 
s  s 


P2X+2  fl2X+l 

2  r  s  cos  n6  "  cos  nPs 

'lrJp2X+l  cos  0  -  y/bg  d0  ’ 
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Pn ’tf;  R2X-1>  b*  P„“t8,tf  iR2X-l> 


*(8), 


,2\ 


g**"'  2\ 

^  />  Ks  cob  n0-  cos  nPg 

^Jp2^"1  cos  0  -  y/be  d0. 


Qnl<'’:aX+l,R2Xl  ^  Qn(S>ly;aX+l'R2x) 


=  11 


n2^+^  2X^1 

'  s  (bgCOS  0-a^)(co3  nB-coe  nP0  *  ) 

if  1  ♦  (bsC08  6  -  aul)(ax  j-  y) 


Q!,9l(y:aXiR2X>  ^  Qn<8)<',;aX'R2X> 

,2X+1 


0  2X. 

2  P  B  (bgCOB  0  -  a^)(cos  n0  -  cob  npfl  ) 


*  ‘  Pf^  1  +  (bgcos  0  -  a^Ka^  -  y) 


d0 


together  with  similar  substitutions  in  b^  <  y  <  b^  \  =  1,  2,  .  . , ,  J. 
After  these  substitutions  are  made,  the  downwash  w(y)  then  becomes 
bounded  everywhere  over  the  wing.  There  are  2J  conditions  of  the 
form  (67)  in  total. 


5,  Application  of  Boundary  Condition  on  the  Wing 

The  same  assumptions  for  the  effect  of  slipstream  rotation  as 
stated  previously  in  Section  B4  will  be  retained  for  the  present  case  of 
multiple  jets.  Since  w(y)  is  even  in  y,  the  boundary  condition  need  be 
applied  for  positive  y  only.  Hence  if  a(y)  is  the  local  geometrical 
angle  of  attack  at  y  of  the  3/4  chord  line,  then,  by  the  Weissinger 
method,  the  required  boundary  condition  becomes 


d0  , 
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wQ(y)/Vj  =  |JLQ(y)  for  b2x+i<  y  ^  b2\»  x  =  o,  1,  .  .  . .  J  ,  (68a) 

(y-a  My-aJ 

wj,  \<yl/Vj  =  Q(y)  4 - Vj -  for  b2X<  y  S  b2X-l  • 

X  =  1,  3,  .  .  . ,  J.  (68b) 

where  to  is  an  even  function  of  its  argument  and  represents  the  slipstream 
rotation,  the  convention  for  the  sign  of  to  being  already  specified  in  Section  B4. 

2J  /  g  i 

Here  we  have  Z  N  unknown  coefficients  A'-,  '  , ,  which  must 
n  s  2n+l 

s  =0 

also  obey  the  2J  constraints  (65)  and  2J  constraints  (67).  Therefore  we 
may  choose  appropriately  (ZNg-  4J)  points  on  the  3/4  chord  line  at  which 
we  apply  condition  (68). 
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III.  COMPUTATIONAL  RESULTS 


The  problem  of  verifying  the  theory  is  two-fold.  First,  an  in¬ 
volved  procedure  for  electronic  machine  computations  must  be  developed 
to  produce  numerical  results  to  apply  the  theory.  Second,  experimental 
data  providing  spanwise  lift  distributions  for  applicable  wing  configurations 
is  required. 

Major  strides  were  taken  in  the  computational  programming  of 
the  theoretical  multiple  jet  relations  developed  in  this  paper.  This  pro¬ 
gram  was  designed  for  the  IBM  Type  709  electronic  computer,  and 
applied  at  the  UCLA  facility.  It  should  be  recognized  that  the  mathematical 
formulation  and  the  computing  procedure  are  highly  interrelated.  In  de¬ 
veloping  the  mathematical  theory,  the  approach  and  selection  of  system 
coordinates  must  be  compatible  with  the  available  computing  techniques. 
However,  all  of  the  computational  difficulties  cannot  initially  be  foreseen 
and  this  results  in  various  modifications  in  the  mathematical  form  of 
presentation  during  the  development  of  a  satisfactory  computational 
procedure. 

The  mathematical  formulation  presented  in  this  paper  yielded 
slight  numerical  inconsistencies  on  several  occasions.  It  is  felt,  how¬ 
ever,  that  the  trends  exhibited  by  the  computed  results  are  reasonable. 

An  improved  formulation  which  shows  promise  of  relieving  some  of  the 
numerical  difficulties  is  presently  under  investigation.  Much  of  the  com¬ 
putational  program  utilized  in  this  paper  would  be  applicable  to  the  new 
formulation. 

Suitable  experimental  data  for  verification  of  the  theory  is  sadly 
lacking.  In  September  of  1959,  NASA  (Ames  Research  Center,  Moffett 
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Field,  Calif.  )  published  an  account  of  force  measurements  of  a  multi-jet 
configuration^*^.  Limited  spanwise  measurements  were  made  during 
the  test,  the  data  for  which  was  not  published  in  their  report.  These 
Ames  spanwise  measurements  were  obtained  upon  request  but  were  not 
suitable  for  correlation  with  the  theory.  Nevertheless,  the  basic  geometry 
of  their  model  was  employed  in  an  example  computation  based  on  the 
present  theory. 

The  Ames  test  was  not  primarily  concerned  with  spanwise  loading 
and  thus  only  a  limited  number  of  spanwise  measurements  were  taken. 

The  propeller  slipstream  velocity  over  the  wing  was  also  not  measured. 

The  effect  of  the  fuselage,  engine  nacelles,  and  asymmetric  loading 
between  the  right  and  left  wing  panels  (spanwise  measurements  were  made 
on  the  left  wing  panel  only)  also  invalidate  any  spanwise  comparison  with 
the  theory.  Also  integration  of  the  Ames  spanwise  lift  data  did  not  agree 
with  their  total  lift  measurements  for  reasons  enumerated  above. 

Fig.  4  presents  theoretical  curves  for  the  spanwise  lift  distribution 
of  the  selected  wing-jet  arrangement  (similar  to  Ames  configuration).  The 
aspect  ratio  is  10  and  the  taper  ratio  is  0,  5.  The  angle  of  attack  is  6°  at 
the  root  chord  and  0°  at  the  tip.  The  spanwise  location  of  the  jets  (pro¬ 
pellers)  is  shown  on  the  sketch  at  the  top  of  the  figure.  The  contribution 
of  the  Bessel  function  terms  for  this  computation  was  found  to  be  small  and 
was  accordingly  disregarded. 

The  most  significant  trend  observed  in  the  theoretical  curves  is 
the  increase  in  the  local  lift  coefficient  near  the  wing  mid- span.  This 
result  is  of  primary  importance  for  structural  wing  design. 

The  curves  also  show  that  lift  is  carried  across  the  narrow  strip 
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between  closely  spaced  jets.  Also,  when  the  jets  are  close  together, 
their  effects  combine  to  produce  a  single  peak.  Obviously,  if  the  jets 
are  a  sizeable  distance  apart,  two  distinct  peaks  will  appear. 


The  maximum  local  lift  coefficient  occurs  near  the  jet  closest 
to  the  mid-span.  This  is  reasonable  since  the  local  lift  of  the  wing 
without  jets  increases  toward  the  center.  The  washout  twist  distribution 
also  amplifies  this  effect.  The  curves  also  show  a  significantly  higher 
local  lift  in  the  region  of  the  outer  jet,  particularly  for  the  higher  jet 
strengths  (lower  p  values). 

Significantly,  there  is  only  one  peak  in  each  of  the  curves  and 
this  peak  shifts  towards  the  wing  mid-span  for  high  p  values.  This  is 
explained  in  terms  of  the  relative  strength  of  the  wing  and  jet,  i.  e.  ,  the 
wing  alone  has  its  maximum  local  lift  (peak)  at  the  mid-span  while  the 
jet  creates  a  peak  at  the  jet  center.  Thus,  the  jet  tends  to  pull  the  peak 
toward  its  center  while  the  wing  tends  to  pull  the  peak  toward  the  wing 
centerline,  resulting  in  a  peak  in  between.  Logically,  this  peak  shifts 
toward  the  wing  centerline  for  weak  jets  (high  p  values). 
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IV.  CONCLUSIONS 


Lifting  surface  theory  has  been  generalized  to  establish  a  foundation 
from  which  many  problems  concerned  with  the  non-uniform  aerodynamics 
of  wing-jet  interaction  may  be  solved.  This  generalization  has  been  applied 
to  the  case  of  a  wing  extending  through  multiple  jets. 

The  multiple  jet  formulation  was  programmed  on  the  IBM  709 
electronic  computer.  Due  to  some  slight  numerical  inconsistencies  an  improved 
formulation  is  presently  under  investigation.  However,  the  trends  indi¬ 
cated  by  the  present  formulation  are  considered  to  be  valid. 

The  spanwise  lift  distribution  of  a  wing  extending  through  four  (4) 
jets  was  computed  (Fig.  4),  A  large  increase  in  lift  inboard  and  within 
the  jets  (propellers)  was  obtained  with  strong  jets  (low  p  values).  Sig¬ 
nificant  increases  in  lift  outboard  of  the  jets  are  also  apparent. 

These  results  provide  a  basis  for  optimizing  multi-jet  wing- 
propeller  configurations,  for  example  by  determining  the  lift  distribution 
of  a  series  of  wings  extending  through  a  systematic  family  of  multiple  jet 
arrangements.  As  demonstrated  in  Ref,  2  the  optimum  planform  for  a 
wing  extending  through  a  single  jet  is  finite  and  :  similar  result  is  antici¬ 
pated  for  multiple  jets. 

The  large  lift  magnification  noted  above  resulting  from  the  wing- 
propeller  interaction  has  wide  applicability  in  the  area  of  V/STOL  aircraft, 
particularly  ir:  the  potential  of  large  improvements  in  STOL  capabilities. 
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V.  RECOMMENDATIONS 


Analytical  solutions  to  the  remaining  two  phases  of  the  four  phase 
fundamental  program  outlined  in  the  Foreword  are  under  way  at  the  present 
time.  With  the  generalization  of  the  lifting  surface  theory  presented  in 
this  report,  the  analytical  solution  of  many  secondary  problems  associated 
with  the  aerodynamics  of  wing-propeller  interaction  can  readily  be  solved. 

It  is  recommended  that  the  analytical  program  be  expanded  to 
encompass  further  secondary  factors  so  that  Rethorst's  basic  solution 
and  the  generalization  of  this  lifting  surface  theory  presented  in  this 
report  are  fully  exploited. 

The  most  pressing  immediate  concern  is  the  need  for  suitable 
multi-jet  experimental  data  for  correlation  with  the  theory,  and  the  need 
for  a  computational  program  of  sufficient  scope  to  delineate  optimum 
configurations  and  gain  an  insight  into  all  the  ramifications  of  the  analyses. 
It  is  recommended  that  an  intensive  experimental  investigation 
of  spanwise  lift  distributions  for  wings  extending  through  multiple  jets 
be  conducted.  The  experimental  program  should  encompass  all  the 
configurations  such  as  the  highly  cambered  wing,  tilt  wing,  etc,,  con¬ 
sidered  in  the  complete  analytical  program  as  outlined  in  the  Foreword, 

It  is  also  recommended  that  the  electronic  computing  program  be 
greatly  expanded  to  fully  utilize  the  analytical  results.  It  is  emphasized 
that  the  computational  program  must  proceed  concurrently  with  the 
analytical  development  since  the  mathematical  formulation  must  be  com¬ 
patible  with  computational  programming  techniques. 
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DEFINITION  OF  ANGULAR  COORDINATES  AND 


FIG. 


FLOW  REGIONS  OF  A  WING  EXTENDING 


SPAN  WISE  STATION 
(NON -  DIMENSIONAL ) 


LOCAL  LIFT  COEFFICIENT, 


SPANWISE  LIFT  DISTRIBUTION 


FRONT  VIEW 


NOTE 

1)  LOCAL  LIFT  COEFFICIENT  BASED  ON 
LOCAL  CHORD  AND  FREE  STREAM 
VELOCITY  OUTSIDE  OF  JET  i_ 

2)  THEORY  BASED  ON  UNIFORM  JET. 
REFINEMENTS  TO  ACCOUNT  FOR 
SLIPSTREAM  ROTATION  ARE  NOT 
INCLUDED  IN  THIS  EXAMPLE  CALCULATION 

3)  THE  SMALL  CONTRIBUTION  OF  BESSEL 
FUNCTION  TERMS  IS  NEGLECTED 

4)  THE  WING  ANGLE  OF  ATTACK  IS  6° 

AT  THE  ROOT  CHORD  AND  O*  AT  THE 
TIP.  THE  ASPECT  RATIO  IS  10  AND  Th'E 
TAPER  RATIO  IS  0.5.  , 


^2..  FREE  STREAM  VELOCITY 


SPANWISE  LOCATION, 


JET  RADIUS 


FIG.  4 
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